Complex Dynamics of a Simple Model for Biological Cells 
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The influence of the internal structure of a biological cell (e.g., a leukocyte) on its dynamics and 
rheology is not yet fully understood. By using a bilamellar vesicle (BLV) consisting of two vesicles 
as a cell model, we find that increasing the size of the inner vesicle (mimicking the nucleus) triggers 
a tank-treading-to-tumbling transition. A new dynamical state is observed, the undulating motion: 
the BLV inclination with respect to the imposed flow oscillates while the outer vesicle develops ro- 
tating lobes. The BLV exhibits a non-Newtonian behavior with a time-dependant apparent viscosity 
during its unsteady motion. Depending on its inclination and on its inner vesicle dynamical state, 
the BLV behaves like a solid or a fluid. These observations provide insight into how a leukocyte 
reacts and adapts its mechanical properties as a response to the flow in microcirculation. 

PACS numbers: 87.16.D-, 87.15. A-, 83.50.Ax, 47.50.-d 



The physical problem — Unilamellar vesicles (ULVs) 
consisting of a single closed phospholipid membrane were 
extensively used as biomimetic model for erythrocytes 
(red blood cells) in the past. They succeeded to repro- 
duce many features known for erythrocytes, like their 
steady shapes in Poiseuille flow [l|, Q or their dynam- 
ical states under shear flow However, for leuko- 
cytes (white blood cells), despite their relevant role in 
the immune system, their dynamics and rheology are 
still poorly understood. The presence of a nucleus al- 
ters their mechanical properties in a non straightforward 
manner In this letter we use a bilamellar vesicle 
(BLV) as a model for biological cells, in particular, a 
leukocyte. While a ULV consists of a single vesicle, a 
BLV consists of two vesicles: an outer larger vesicle (the 
cell) enclosing an inner smaller one (mimicking the nu- 
cleus) , see Fig. [TJi. For the first time, we study nu- 
merically the dynamics of a BLV under shear flow and 
we investigate how the dynamical and rhcological prop- 
erties of a leukocyte (modeled by a BLV) are affected by 
varying the size and the deformability of the nucleus as 
well as the amount of fluid enclosed between the nucleus 
and the cell. We demonstrate that leukocytes cannot be 
approximated simply with a fluid-filled particle enclosing 
a homogeneous fluid without an internal structure as it 
has been used, for example, in Refs. [1,0]. The reason 
is that leukocytes adapt their mechanical properties and 
act as a solid or as a liquid depending on how they are 
deformed by the imposed fluid [(§ . 

Computational method — We consider two concentric 
vesicles in two dimensions and place them in a shear flow 
generated between two parallel plates. We designate by 
i?out and R m < P ou t the effective radii of the outer and 
the inner vesicle (in 2D, R = P/(2tt), where P is the 
vesicle perimeter) . All fluids are considered to be incom- 



pressible, Newtonian and of the same viscosity rj. Their 
flow is solved by the lattice-Boltzmann method and the 
fluid-vesicle two-way coupling is achieved employing the 
immersed boundary method (see for details). Both 
vesicle membranes are locally inextensible and experience 
resistance to bending with the same rigidity k. They ex- 
ert a reaction force 
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f= K^ + 2j- aT T + ^ t ' (1) 

on the surrounding fluid where n and t are the unit nor- 
mal and tangent vectors, c is the local curvature, s is the 
curvilinear coordinate and a* is the local effective sur- 
face tension. We do not use additional interaction forces 
between the membranes; they interact only hydrodynam- 
ically. The distance between the plates is chosen as such 
that the effect of wall confinement is negligible El ■ 
First, we investigate how the dynamics of the BLV is af- 
fected by varying two parameters: i) Pi n , to study the 
effect of its internal structure and ii) the capillary num- 
ber Ca = rjjR^ ut /K, where 7 is the shear rate, to study 
the effect of the membrane deformability. Second, we 
investigate how the dynamics of the inner vesicle is af- 
fected by the flow induced by the outer vesicle by varying 
the swelling degree A ou t defined as 4tt A out /P 2 ut (A>ut is 
the outer vesicle area) while keeping all other parame- 
ters fixed. All simulations are performed in the Stokes 
regime: O(Rc) = 1CT 2 , where Re = pjRl ut /n is the 
Reynolds number and p is the fluid density. 

Observations and discussion — Both vesicles are de- 
flated and have the same swelling degree, A; n = A out = 
0.9. We vary P; n while keeping P ou t fixed. The ob- 
tained dynamics is shown in Fig. Q] as a function of 
< Pin/Pout < 1- For smaller inner vesicles (Fig. QJa)- 
(c)), both vesicles perform a steady tank-treading motion 
(they assume a steady inclination angle with respect to 
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FIG. 1: Snapshots showing the dynamics of a bilamellar vesicle (BLV) in shear flow at low deformation regime (Ca = 0.5). 
The ratio Rin/Rout is given below each image. For small inner vesicles (a)-(c), the BLV performs steady tank- treading. The 
angle between the main axis of the BLV and the flow direction 9 out decreases with R[ n (see Fig. [2}. Beyond a threshold of Ri n , 
the BLV starts tumbling (d). This is similar to the tank-treading-to-tumbling transition observed for viscous vesicles [EL ITll| . 
but is triggered here (in the absence of a viscosity contrast) by the presence of the inner vesicle. 



the flow while their membranes undergo a tank-treading 
like motion). The BLV aligns more and more with the 
flow when i?i n is increased. Due to the symmetry, both 
vesicle centers are stationary. Beyond a threshold of 
Ri n , the motion of the BLV transits from steady tank- 
treading to unsteady tumbling motion (rotation as solid 
elongated particle). Fig.QJd) shows snapshots of a tum- 
bling BLV with i?in/-Rout = 0.75. The inner vesicle as- 
sumes a relative angle with respect to the main axis of 
the outer one. During the tumbling motion, the mean 
value of # ut (the angle defined by the main long axis of 
the BLV and the flow direction) is zero. When plotting 
#out as function of Rin/Rout (Fig. 0, we see that 8 out 
decreases with increasing Ri n /R on t until it vanishes at 
a critical value. Then, tank-treading (TT) is replaced 
by tumbling (TB). The TT-TB transition is known for 
viscous ULVs where it is induced by increasing the vis- 
cosity contrast A (ratio between the internal and external 
fluid viscosities) beyond a given threshold [ll| . Here, 
however, all fluids have the same viscosity and the tran- 
sition is induced solely by the presence of the encap- 
sulated vesicle and by enlarging its size. Veerapaneni 
et al. [l2l ] predicted a similar transition for non-viscous 
ULVs enclosing a solid particle. They claim that the 
inclusion increases the apparent internal viscosity lead- 
ing to the transition as observed for inclusion-free ULVs 
with A > 1. To investigate the effect of the apparent 
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FIG. 2: The inclination angle 9 out (left axis) and the apparent 
viscosity contrast A* (right axis) as a function of the radius 
ratio. An increase of Rin/Rout induces an increase in A* and 
the BLV becomes more and more viscous. Consequently, 9 out 
decreases until it vanishes at the transition point. 
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FIG. 3: The inclination angle vs. the viscosity contrast for a 
viscous ULV and a BLV: the angle decreases with increasing 
A or A*. However, for a given A or A*, the angle of the 
BLV is larger than that of the viscous ULV demonstrating a 
difference in the quantitative behavior of both. 



internal viscosity, we compute rj* = 77 + (<J xy ) / (S X y) , 



where (a xy ) = — 



TdSli 



ds (fxr y )/Ai n is the average ex- 
cess shear stress in the inner vesicle domain fl m and 
(Sxy) = § dn ds (n x u y +n y u x )/A out is the average shear 
rate in the outer vesicle domain f2 ou t. Here, r and u are 
the position and the velocity of a membrane point, re- 
spectively. The apparent viscosity contrast of the BLV 
A* = V*/ 1 ! vs - -Rin/-Rout is depicted in Fig. [2] increasing 
Rin leads to a monotonic increase of A*. The TT-TB 
transition for the BLV takes place at a critical value of 
A* r = 6.9. This value is close to the critical viscosity 
contrast A cr = 6.6 required for a viscous ULV to un- 
dergo the same transition (for the same swelling degree 
of 0.9). A systematic comparison of the inclination angle 
(6, Gout) vs. the viscosity contrast (A, A*) of a viscous 
ULV and a BLV (see Fig. [3]) shows that both exhibit 
similar qualitative behavior: the angle decreases with in- 
creasing viscosity contrast until it vanishes at the transi- 
tion point. However, for all viscosity contrasts, the angle 
of the BLV is found to be larger than that of the viscous 
ULV (6 ou t > 6), especially at larger A* corresponding to 
larger inner vesicles. This demonstrates that a BLV does 
not behave exactly as a viscous ULV for which the in- 
ternal fluid is a homogeneous medium. Moreover, when 
tumbling, the apparent internal viscosity 77* of the BLV 
is a time-dependent quantity and varies in a coherent 
manner with 9 ou t, see Fig. 2J It diverges to the limit of 
a solid medium when 6 ou t = ±7r/4 (direction of the elon- 



3 



gation/comprcssion of the shear flow) since the average 
shear rate (S xy ) vanishes at that point. This is a sig- 
nature of the non-Newtonian rheological behavior. The 
BLV internal medium changes its apparent viscosity as a 
response to the orientation with respect to the flow. 
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FIG. 4: Time evolution of the inclination angle 8 ou t (in de- 
grees) and the apparent internal viscosity rf of a tumbling 
BLV (Rin/Rout = 0.75, Ca = 0.5). The data corresponds to 
the snapshots in Fig. 1(d). The time evolution of the shear 
stress {(Txy) and shear rate (S X y) is also shown (in numerical 
units). The time dependence of n* suggests non-Newtonian 
fluid properties. Note especially rf — ► oo (solid limit) when 
6 = ±vr/4 (i.e. (S xv ) -> 0). 



At higher Ca the membrane deformability becomes im- 
portant and thus the BLV deforms substantially. For a 
smaller inner vesicle, the BLV again tank-treads. How- 
ever, for larger R ln , in contrast to the limit of small de- 
formation, we surprisingly observe that the BLV does not 
tumble anymore. It rather performs a new type of un- 
steady motion (Fig. [5]). The inner vesicle undergoes a 
swinging motion; the main axis oscillates about a posi- 
tive mean angle (#;„ in Fig.H^a)) while its shape does not 
deform. Such motion is known for red blood cells and 
capsules The outer vesicle exhibits a non-regular 

motion: although its main long axis performs oscillations 
about a positive mean angle as well ((9 out in Fig. [5ja)), 
its shape undergoes larger undulations, i.e., its membrane 
buckles. It develops two oscillating lobes for intermedi- 
ate sized inner vesicles or four rotating lobes for larger 
inner vesicles (Fig.JSJb)). Both unsteady motions cannot 
be qualified as vacillating-breathing (VB) [Til ], swinging 
(SW) or trembling (TR) [l5|. An almost similar fea- 
ture as shown in Fig. [5] has been recently observed ex- 
perimentally by Pommella et al. [l6| for a surfactant 
multilamellar droplet subjected to strong shear. The au- 
thors describe the droplet motion as VB, but we disagree 
with this classification since the angle is found to oscil- 
late about a non-zero mean value as it does for SW. Yet, 
SW can be ruled out as well since it is observed in the 
small deformation limit. The dynamic mode of the BLV 
can neither be described as TR which indeed is charac- 




FIG. 5: (a) Oscillations of the main axes of the inner (0i n ) 
and outer (8 out ) vesicles of a BLV during the undulating mo- 
tion (Ri n /Rout = 0.85, Ca = 10). (b) Snapshots, taken at 
equal time intervals, showing four lobes of the outer vesicle 
membrane and their rotation. The straight line denotes the 
consequent locations of two opposite lobes. The undulating 
regime replaces tumbling for larger i?i n /iiout and larger Ca. 



terized by the formation of lobes [15( but also requires 
the shape to become perfectly elliptical at a certain mo- 
ment. This is impossible for BLV due to the presence of 
the large inner obstacle. The appearance of this new un- 
steady motion (for larger Ca and larger Ri n /R ou t) that 
we name undulating motion cannot be explained solely 
based on the apparent viscosity contrast argument. The 
inner vesicle disturbs the motion of the outer one. By 
increasing its size, the thickness of the fluid layer be- 
tween the membranes decreases to become a thin liquid 
film. The outer membrane tries to tank-tread under the 
external applied shear. However, the presence of the in- 
ner vesicle prevents this and thus it slides over the inner 
membrane, which plays the role of a nearly solid obsta- 
cle. Its shape is less deformablc because i) the inner vesi- 
cle is smaller Ca oc R 3 , and ii) the outer vesicle shields 
the inner one from the external flow. A thorough un- 
derstanding of the appearance of the undulating motion 
is still missing. However, a relation to the Marangoni 
effect can be proposed: at large deformations, bending 
becomes less important than tension. We observe a non- 
uniform distribution of the surface tension for the BLV, 
da* /ds 7^ 0. This is in line with the surfactant multil- 
amellar droplet (l6j and the instability of thin liquid films 
on a solid substrate 



So far we described how the inner vesicle alters the 
dynamics of the outer vesicle and so of the BLV. Next, 
we examine how the outer vesicle in turn influences the 
dynamics of the inner vesicle. This is controlled by the 
amount of the fluid between the two membranes. To un- 
derstand its relevance we vary its amount by swelling 
(adding fluid) or deflating (removing fluid) the outer 
vesicle. We vary A out between 0.9 and 1, keeping all 
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(a) A out = 0.90 (b) A out = 0.98 (c) A out = 1.00 




(d) T/S = 11.56 (e) T/S = 27.62 (f) T/S = 243.90 



FIG. 6: The dynamical states of the inner vesicle as a func- 
tion of the outer vesicle swelling degree A out . (a) Tank- 
treading motion, (b) tumbling while the outer vesicle performs 
a breathing-like motion, and (c) tumbling with a constant an- 
gular velocity. Panels (d)-(f) show the resulting flow field (in 
absence of an inner vesicle). It evolves from mixed flow to 
pure rotational flow when A out — > 1. T/S denotes the ratio 
of flow vorticity and shear rate. 



other parameters unchanged. Consequently, we observe 
another TT-TB transition, but this time for the inner 
vesicle as shown in Fig. [BJa)-(c). For A = 0.9, both 
vesicles tank-tread (Fig. Eta)). Above a critical value of 
A ou t = 0.98, the inner vesicle starts tumbling while the 
outer experiences a breathing- like motion (Fig. Ob) ) . For 
A out = 1 (Fig. Etc)), the BLV behaves exactly as a solid 
body. The two vesicles rotate with the same angular ve- 
locity and also the fluids, between the membranes and 
within the inner membrane, behave like a solid medium. 
The dynamical transition observed here for the inner vesi- 
cle is expected to modify the apparent viscosity within 
the BLV. A similar link between the rheology and the 
micro-dynamics has been observed for red blood cells 
[l8j |. Although the BLV dynamics for the three cases 
(Fig. Ela)-(c)) is apparently similar for an outside ob- 
server, their rheological properties may differ due to the 
dynamical state (TT or TB) of the inner vesicle. The 
dynamical transition induced by varying A ou t alone can 
be explained by the theory of Lebedev et al. who 
predicted that even a no n- viscous vesicle (A = 1) can 
undergo a TT-TB transition by increasing the rotational 
component of the external imposed flow. This was later 
confirmed experimentally by Deschamps et al. pol ]. In 
our case, the inner vesicle is subjected to the flow in- 
duced by the tank-treading motion of the outer mem- 
brane. Fig. Eld)-(f) depict the generated undisturbed 
flow (in the absence of the inner vesicle) for each A ou t- 
A pure rotational flow is obtained for A ou t = 1- For 
A ut = 0.9 and A ou t = 0.98, a mixed flow is gener- 
ated, i.e., a combination of pure shear and pure rotational 



flows. We quantify the relative importance of the rota- 
tional and the elongational components using the quan- 
tity T/S, where T = | <f dn ds(n x u y — n y u x ) is the 

vorticity and S = ^\/^{(S X y) 2 ) is the shear magnitude. 
T/S = 1 for pure shear flow and T/S — > oo for pure ro- 
tational flow. The computed values T/S for each A ou t 
are shown in Fig. |B{d)-(f). The TT-TB transition is 
induced by increasing T/S beyond the critical value of 
27.62 (corresponding to the critical value A ou t = 0.98). 
Using Lebedev's parameters, S = 77rCa/ 'y/36 and A = 
4(1 + 23A/32)V^/V307r (5 is the excess perimeter in 2D 
or excess area in 3D), the transition occurs in our case 
at S CI = 0.42 and A cr = 11.41. This A cr is larger than 
expected when compared to the case of a viscous vesicle 
under shear flow (A cr = 1.38) II ill or a non- viscous vesi- 
cle in general flow (A cr = 1.2) [20(. At low deformation 
(S -> 0), A cr is_independent of S anyway [H Ho[. The 
discrepancy in A cr cc A can be explained by taking into 
account that Lebedev's theory has been formulated for 
unbounded flows while, in our case, the inner vesicle is 
strongly confined by the outer vesicle membrane. Con- 
finement is found to shift the critical point of the TT-TB 
transition to larger values of A [llj. Thus, Lebedev's the- 
ory can explain only qualitatively but quantitatively the 
TT-TB transition observed for the inner vesicle in Fig. [6] 

Summary & conclusions — We showed that even a 
non-viscous vesicle (known to perform only TT under 
shear flow) exhibits rich complex dynamics when it en- 
capsulates another non-viscous vesicle. Increasing the 
size of the inner vesicle triggers a dynamical transition 
from TT to TB or the newly found motion: undulat- 
ing. The BLV internal medium (the fluid plus the inner 
vesicle) displays non-Newtonian behavior with a time- 
dependent apparent viscosity during unsteady motions: 
the same BLV behaves like a solid or a fluid depending 
on its orientation and the dynamical state of its inner 
vesicle. Our results suggest that a leukocyte cannot be 
simply mimicked with a solid spherical particle or with 
an inclusion-free vesicle enclosing a homogenous Newto- 
nian fluid. The presence of an internal structure dictates 
its dynamical and rheological response to imposed flow. 
Further investigation of the consequences of our results 
on the margination and adhesion of leukocytes in the 
microcirculation, and its physiological and pathological 
implications, is to be considered in the future. 
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